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Abstract 

An experimental study of the production of up-going charged particles in inelas- 
tic interactions of down-going underground muons is reported, using data obtained 
from the MACRO detector at the Gran Sasso Laboratory. In a sample of 12.2 x 10 6 
single muons, corresponding to a detector livetime of 1.55 y, 243 events are ob- 
served having an up-going particle associated with a down-going muon. These 
events are analysed to determine the range and emission angle distributions of the 
up-going particle, corrected for detection and reconstruction efficiency. 

Measurements of the muon neutrino flux by underground detectors are often 
based on the observation of through-going and stopping muons produced in i/„ 
interactions in the rock below the detector. Up-going particles produced by an 
undetected down-going muon are a potential background source in these measure- 
ments. The implications of this background for neutrino studies using MACRO 
are discussed. 



1. Introduction 



One of the primary goals of the MACRO detector at the Gran Sasso underground labo- 
ratory is the measurement of the flux of atmospheric muon neutrinos. In this paper we 
present the measurement of a potential background in these studies; charged up-going 
particles (primarily pions) produced in inelastic interactions of down-going muons in 
the rock surrounding the detector. If the down-going muon is not detected the up-going 
pion may be misidentified as resulting from a neutrino-induced event. This will occur, 
for example, when the n passes near but not through the detector. This type of event, 
therefore, represents a background in the measurement of the flux of neutrino-induced 
upward through-going and stopping muons. The events under study, although similar 
in topology to muon backscattering, have an observed rate which is many orders of 
magnitude larger than that expected from the large angle scattering of muons [I], |2|]. 
The relatively low energy up-going tracks are primarily charged pions produced in hard 
muon scattering on a nucleon N, /i + iV — > /j, + n + X. Hadron production by un- 
derground muons has been recently discussed || in conjunction with the problem of 
neutron background in the measurement of the v e /^^ atmospheric neutrinos ratio by 
large water Cherenkov detectors M\ . Hadron production by muon interactions in matter 
is also a source of background in radiochemical solar neutrino experiments || . 

In this paper we present the results of a study of these events carried out with the 
MACRO detector. In Section 2 the MACRO apparatus and the method used to identify 
muon-induced up-going charged particles are described. The details of the analysis of 
the data are given in Section 3. In Section 4 the observed range and emission angle 
distributions for the up-going particle are presented. In Section 5 the detector response 
and the tracking efficiency for these events are calculated. In Section 6 the results are 
interpreted within the framework of a model of the photonuclear interactions of high 
energy muons provided by the FLUKA generator. Finally, in Section 7 the background 
due to up-going particles produced by hard muon scattering in the measurement of 
neutrino- induced muons is calculated for the MACRO experiment. 

2. Muon-induced Up-going Particles in MACRO. 

The MACRO detector, described in detail in Ref. |J, is located in Hall B of the Gran 
Sasso underground laboratory. The detector has global dimensions of 76.6m x 12m x 
9.3m, and is divided longitudinally in six similar supermodules and vertically in a lower 
(4.8 m high) and an upper part (4.5 m high). The active detectors include horizontal 
and vertical planes of limited streamer tubes for particle tracking, and liquid scintillation 
counters for fast timing. In the lower part, the eight inner planes of limited streamer 
tubes are separated by passive absorber (iron and rock ~ 60g cm~ 2 ) in order to set a 
minimum threshold of ~ 1 GeV for vertical muons crossing the detector. The upper 
part of the detector is an open volume containing electronics and other equipment. The 
horizontal streamer tube planes are instrumented with external pick-up strips at an 
angle of 26.5° with respect to the streamer tube's wire, providing stereo readout of the 
detector hits. The transit time of particles through the detector is obtained from the 
scintillation counters by measuring the mean time at which signals are observed at the 
two ends of each counter, and taking the difference in the measured mean time between 




RUN = 9967 EVENT= 3941 1 1-MAY-95 07:36:38 

Figure 1 : On-line display of a typical muon interaction in the rock below the apparatus giving 
rise to an up-going charged particle. The rectangular boxes indicate scintillator counter hits 
and the points are streamer tube hits. The tracks are shown in the wire view of the streamer 
tubes. Only two supermodules are drawn. 

upper and lower counters. The particle direction in the detector is given by the sign of 
this difference. 

The apparatus can observe charged particles produced at large angles by hard muon 
interactions if the reaction happens either inside the apparatus or in the rock and con- 
crete below the apparatus. In the latter secondary particle produced at a large 
angle in the reaction is seen as an up-going particle in the apparatus (see Fig. [I]). If 
the muon passes through the detector, such an event consists of two tracks converging 
somewhere below the detector, both in the wire and strip projective views of the track- 
ing system. The secondary particle reaches the scintillator counters with a positive time 
delay with respect to the counters hit by the muon. 

In the sample of muon plus up-going pion events, the down-going particle is usually 
the longer of the two tracks, and is reconstructed by the standard MACRO tracking 
procedure ||. This tracking algorithm is optimized for the reconstruction of single and 
multiple muon events, and requires that at least four central streamer tube planes are 
crossed by the particle, corresponding to a minimum of ~ 200 g cm~ 2 of detector gram- 
mage traversed. For this study, a tracking algorithm was developed to reconstruct the 
short, diverging tracks characteristic of the events of interest. This algorithm searches 
for alignments between a cluster of scintillation counters and at least two streamer tube 
hits in a 4 meter wide region centered on the counter cluster. The scintillator cluster is 
defined as a contiguous group of neighboring counters simultaneously hit by one or more 
particles. The reconstructed coordinates of the scintillation cluster, obtained from the 
identity of the hit counters and the timing difference between the two ends of the coun- 
ters, are used by the tracking algorithm as an additional track point. A track candidate 
consists of an alignment between at least three hits in a 4 meter wide region, and is 
referred to below as a short track. The minimum vertical depth for events reconstructed 
by this algorithm is approximately 25 g cmT 2 . 




Figure 2: Sketch of an event containing a down-going muon with direction (#„, </?„) producing 
an up- going pion with direction (-3^,^-,^). The arrows indicate the direction of the particles, 
deduced from the timing information obtained from the scintillation counters (Ti,T2 and T3). 
The filled points correspond to the streamer tube hits. The zenith and azimuth emission angles 
Aa, Aip of the pion with respect to the muon direction are indicated. 

3. Event Selection 

The data sample used in this study includes events in which a single muon is recon- 
structed by the standard tracking algorithm (see Fig. |2|, for a sketch). The short track 
algorithm is then used to search for events in which a secondary up-going pion enters the 
detector from below and crosses the bottom layer of scintillator counters and at least two 
streamer tube planes before coming to rest in the detector. Finally, the timing informa- 
tion from the scintillation counters is used to identify the direction of the two particles. 
As demonstrated below, the pion and the muon are resolved with high efficiency if the 
wire hits belonging to the two tracks are separated by at least ~ 75 cm. 

Two different data periods have been analyzed. In the first (sample A, from December 
1992 through June 1993), the upper part of the apparatus was not in acquisition. In 
the second period (sample B, from April 1994 until January 1996) the entire apparatus, 
including the upper part was in acquisition. Table 1 contains the detector livetime, the 
number of single downgoing muons, the number of reconstructed short tracks, and the 
number of final muon plus up-going pion events obtained in the two data taking periods. 
As shown in the table, a short track is reconstructed in 31717 of the 12.2 x 10 6 single 
muon events. These are referred to below as double track events. A visual scan of a 
subsample of 15% of these events has been performed to determine their nature. In the 
most common case (76% of events in sample A, and 69% in sample B), the short track 
in the event belongs to a second downgoing muon, parallel to the first. The difference in 
the fraction of events of this type between the two data samples is due to the improved 
acceptance for double muon events during the period in which the entire apparatus 
was in acquisition. Other topologies, all producing an incorrectly reconstructed short 
track, are due to electromagnetic showers inside the detector (17% of events in the two 



Sample 


Single 
downgoing muons 


Double 

tracks 


Livetime 
00 


Down muon 
+ up pion 


A 


3.3- 10 6 


12216 


4235 


65 


B 


8.9- 10 6 


19501 


9351 


178 


A+ B 


12.2 • 10 6 


31717 


13586 


243 



Table 1: Number of single down-going muons and of events with an additional short track 
(double tracks) in the two data periods A and B. In the fourth column the livetime (h) 
of the apparatus is shown. The number of identified muon plus up-going pion events is 
given in the last column. 

periods), incorrect track point assignment (4%), muon interactions occurring inside the 
detector (5%) and low momentum muons undergoing significant multiple scattering in 
the detector (2%). In all of these event types the two reconstructed tracks do not possess 
a common vertex point below the detector. 

The event topologies identified as background in the visual scan are rejected by 
applying the following analysis cuts to the sample of double track events. Double muon 
events are rejected by requiring that the direction of the two tracks differ by more than 
26°. Electromagnetic showers are rejected by requiring that the total number of streamer 
tube hits not belonging to the muon track be less than 20, while muon interactions inside 
the apparatus and low energy muons are rejected by requiring that the intersection point 
of the two tracks be outside the detector volume. 3467/31717 events pass these cuts and 
have been visually scanned to obtain the final data sample. Events eliminated in this 
final visual scan include in most cases multiple and/or wrongly reconstructed short 
tracks, produced by a double muon event or by a muon electromagnetic shower. As a 
result, 65 events in sample A and 178 in sample B are classified as muon plus up-going 
pion events. 

The probability of rejecting a genuine muon plus up-going pion event through the 
software cuts has been evaluated by scanning a random subsample of 6200 double track 
events, representing ~20% of the full double track event set, for muon plus up-going 
pion events. 54 events are found, of which 3 are rejected when the software cuts are 
applied. The corresponding selection efficiency is e s = 94%. An independent check 
of e s is obtained by applying the software cuts to the sample of simulated (discussed 
below) down-going muons plus up-going pions; 151/3347=4.5% of the simulated events 
are discarded. Based on the number of events in the final data sample, and a selection 
efficiency e s = 94%, the rate of detected up-going pions per down-going muon at the 
MACRO depth is: 

243 ... 
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time delay between counters (ns) 

Figure 3: Total path length (see text) versus time delay At = (T 3 — T 2 ) between scintillation 
counter hits in the bottom plane. Only the 168 events for which \z™ — z%\ < 50 cm are included. 
The large spread around the central line is due to the combined uncertainty in the total path 
length and timing reconstruction. As an example, error bars are drawn for two events. For 
three events, T 2 is missing and At = 

4. Analysis 

The 243 events in the final data sample are now used to obtain the distribution of the 
range of the up-going pion, the angle between the /1 and tt, and the radial distance D 
between the pion and the muon at a vertical position corresponding to the center of the 
bottom layer of scintillation counters. 

Referring to Fig. |2|, the point at which the muon interaction occurs^ Xu — (^x u y y^jj z^j^ 
is defined as the intersection of the observed muon and pion tracks. The vertex depth z v 
is obtained separately in the wire (z™) and strip (z*) projection. In 75 of the 243 events 
in the final data set the difference between the vertex depth in the wire and strip views 
exceeds 50 cm; this event subsample is denned as poorly reconstructed in the following. 

The timing information obtained from the scintillation counter system is used to 
determine the direction of motion of each particle, and as an independent check of the 
vertex location obtained from the tracking system. Using the convention of Fig. ^|, 
a down-going particle will have a transit time difference (T 2 — T{) > 0. The bottom 
counter crossed by the short track is reached at a time T 3 , with (T 3 — T 2 ) > 0. For 206 
of the 243 events in the final data sample the requirement (T 3 — T 2 ) > is fulfilled. For 
the remaining 37 events, either the down-going muon misses the counters of the bottom 
scintillation plane (T 2 not present) or the two particles intersect scintillation counters 
in supermodules read out by different data acquisition systems (for a small fraction 
of the time, the electronics to synchronize the information between the supermodule 
readout systems was non-functional). Because the timing error seen in these events can 
be accounted for, all 37 are classified as muon plus up-going pion events and are included 
in the subsequent analysis. 

The total path length between the two scintillation counter hits, given by the path 
length of the muon from the scintillator T 2 to x v plus the path length of the pion from 
x v to the scintillator T 3 , is determined for events in the final data sample which satisfy 
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Figure 4: a) Distribution of the reconstructed total range Rt of the up- going particles for 
which \z™ — z%\ < 50 cm. Events with a poorly reconstructed vertex are excluded from this 
distribution because of the large uncertainty in Rt for these events, b) Distribution of the 
grammage traversed in the detector, Ri, for all the 243 up- going particles in the final data set. 



the condition \z™ — z*\ < 50 cm. The correlation between the total path length and 
the time delay (T 3 — T 2 ) is shown in Fig. |^. The linear correlation seen in this figure 
is consistent with the hypothesis that the up-going particle is produced at the vertex 
point. The uncertainty in the total path length depends on the uncertainty in the vertex 
location, and ranges from 50 cm up to 250 cm for the worst. The average uncertainty 
in the time delay between the two scintillation counters is approximately 2 ns for this 
data set. 

The range of the up-going particle, R t , is denned as the distance between the in- 
teraction point (x v ) and the position of the uppermost streamer tube hit in the track, 
(x s ). The range obtained in this way is, in most cases, an underestimate of the true 
range of the pion, by an amount as large as the vertical depth of the concrete absorber 
layer between streamer tube planes, 60 g cm~ 2 . The particle range is calculated by 
dividing the path length from x v to x s into 200 steps and considering the density of the 
material at each step. The sum of the product of the density times the step length gives 
the total range in g cm~ 2 . This procedure makes use of a detailed description of the 
distribution and density of the material inside and around the apparatus. In Fig. [|a the 
distribution of the total range of the pion, R t , is shown. Only muon plus up-going pion 
events satisfying the condition \z™ — < 50 cm are included in this distribution. The 
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Figure 5: a) Measured distribution of the ir — \x scattering angle Aa and b) of the distance 
D between the down-going muon and the up-going pion at the z location of the center of the 
bottom layer of scintillation counters. 



partial range, R i} of the up-going particle inside the apparatus is also calculated, and is 
shown in Fig. |]b. Ri is calculated from the bottom of the detector to the last active 
streamer tube hit in the up-going track. All the 243 muon plus up-going pion events 
are included in Fig. |]b, since the position of the vertex is not used in the calculation of 
Ri. It can be seen from these figures that a significant fraction of these events contain 
an up-going particle which penetrates enough material to be incorrectly identified as an 
up-going muon when the down-going particle is undetected. 

Fig. [5]a shows the measured distribution of the pion emission angle, Aa, while Fig. 
[5|b shows the distribution of the distance D between the muon and the up-going pion 
at a vertical position corresponding to the center of the bottom layer of the scintillator 
counters. The distributions, as discussed in the next section, decrease for small D as 
a result of a low event reconstruction efficiency, and for small Aa as a consequence of 
the low muon rate at large zenith angle. The radial distance between the muon and the 
up-going pion is less than 4 m, corresponding to roughly one third of the lateral size of 
one MACRO supermodule, in 90% of the events, as discussed in Section 7. 
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Figure 6: Distribution of the probability A for events to be accepted, of the tracking efficiency 
E, and of the reconstruction efficiency AE versus (a) Aa , (b) Aip , (c) D, and (d) R n . Each 
quantity is averaged over the three remaining variables. Full points: A; empty points: E; solid 
curves: AE. 



5. Reconstruction Efficiency 

The probability of detecting an up-going pion produced by a down-going muon depends 
on the emission angles of the pion with respect to the muon, the pion energy, and the 
vertex depth. To determine the detection probability for these events in MACRO, a 
sample of 5.0 x 10 5 simulated muon plus up-going pion events is distributed uniformly 
over cells in the parameter space V defined by the variables Aa, Ap, z v and the pion 
range R v (see Fig.^). The zenith angle Aa of the pion with respect to the incoming 
muon direction is constrained to lie in the range 0° < Aa < 180°. The vertex location 
x v is taken at random locations below the detector, with the vertex depth at a distance 
5 < d v < 205 cm from the bottom layer. The pion range inside the lower detector 
is chosen between < R n < 1000 g cm~ 2 . This corresponds to a value of the pion 
momentum at the detector surface, p n , in the interval 0.1 < p n < 2.0 GeV/c. In each 
simulated event, the down- going muon direction ip^) is extracted from the measured 
distribution of atmospheric single muons. The point at which the muon enters the 
detector is distributed uniformly over the detector surface. The pion direction, ($„., p T ), 
is given by $ n = 180° — $ M — Aa and p T = p^ — Ap. Using the events uniformly 
distributed on respect these parameters, and the actual angular distribution of the muons 



at MACRO depth, the reconstruction efficiency, averaged on the muon directions, has 
been obtained as a function of the chosen parameter set. 

G MACRO §, a GEANT-based § Monte Carlo of the MACRO detector, is used 
to model the response to these particles. It describes the experimental apparatus (i.e. 
geometry, density, and the detector response) in all its details. The propagation of low- 
energy pions is potentially sensitive to the treatment of charge-exchange and in-flight 
absorption reactions. For this reason, the propagation of pions in the detector using the 
GEANT-GHEISHA (G-G) hadronic interface has been compared with results using the 
GEANT-FLUKA (G-F) interface (see 0] for references). We find that the two models 
give essentially the same detection probability for muon plus up-going pion events. For 
example, G-F predicts that a 500 MeV pion entering the detector from the bottom 
layer has a 20% probability of a range within the detector smaller than 25 g cm~ 2 , while 
G-G predicts a probability of 23%. This range corresponds to the minimum distance 
required to identify the up-going pion, and as a consequence, the number of accepted 
events does not differ significantly for the two cases. On the other hand, the energy 
distribution of the up-going pions, unfolded from the measurement of the range, is quite 
different for the two models. This can be illustrated by noting that the same 500 MeV 
pion has a probability of 57% of having a range smaller than 100 g cmT 2 in G-F, while 
G-G predicts a 74% probability. In the following, we report the results obtained with 
the GEANT-FLUKA simulation. 

The result of processing the simulated muon plus up-going pion events in GMACRO 
is an output data set whose format is essentially identical to that of real data. This 
data set has been processed with the analysis chain used for the real data, described in 
section 4. A simulated event is defined as accepted, if the down-going muon is tracked 
and if the pion enters the detector and crosses the lower two horizontal streamer tube 
planes and the horizontal scintillator layer in between. The minimum amount of material 
traversed by a pion satisfying this requirement is ~ 25 g cm~ 2 . The event is defined 
as tracked if the pion track is reconstructed, as described in Section 2. The acceptance 
probability and tracking efficiency for the simulated events is now obtained for each cell 
in the parameter space V = (R n , Aa, Atp, D(z v )). The radial distance D is uniquely 
determined for a given pion and muon directions ($„., (p n ), (i?^, <^ M ) and vertex depth z v . 
The acceptance probability for an event is defined as 

. .„ . . _ number of accepted event , . 

A{R n , Aa, A<p, D) = /. f— 2 

number oj simulated events 

while the tracking efficiency is defined as 

. „ _ number of tracked events , . 

E(R W , Aa, Aip } D) = j — 3 

number of accepted events 

Fig. |6| shows the acceptance probability A, the tracking efficiency E, and the recon- 
struction efficiency AE as a function of Aa , Atp , D and R n . In each case, the values 
of A and E are obtained by averaging over the remaining three variables. Referring to 
Fig. |6], very few pions emitted with a small Aa (Fig. |6]a) are detected, due to the fact 
that these pions must be produced by a near-horizontal muon, and the integrated muon 
flux with "dp > 75° is 6 x 10~ 3 of the total. In Fig. ^|c, the decrease in the tracking 
efficiency for small D is a result of the merging of the pion and the muon tracks. Fi- 
nally, the reconstruction efficiency increases as a function of the pion range within the 
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Figure 7: Reconstruction efficiency AE in the (R n , Aa, D) parameter space, for 4 intervals of 
the distance D, 12 bins of Aa and 10 bins of the measured range of the pion inside the detector 
R n . The computed values of AE are proportional to the box area. For example, the upper 
row in (d) corresponds to values of (AE ■ 100) equal to 1, 1, 3, 18, 29, 48, 52, 47, 40, 30, 48 and 76. 
In each cell the number of detected muon plus up-going pion events at the given (R n , Act, D) 
values is indicated. 



detector (Fig. [5]d), because the tracking efficiency increases with the number of streamer 
tube hits in the track. Since (Fig. ||b) AE is essentially independent of the azimuthal 
scattering angle Aip, it is assumed in the following that AE = AE(R K , Aa, D). 

In Fig. 0, the value of AE is shown as a function of R w and Aa in four intervals of 
the distance D: (a) D < 100 cm, (b) 100 < D < 200 cm, (c) 200 < D < 400 cm and 
(d) D > 400 cm. The box sizes are proportional to the reconstruction efficiency AE 
for detecting an event at a given (R w ,Aa,D). In the same figure, the numbers inside 
the boxes correspond to the distribution of real events. In the four intervals of D there 
are 11, 78, 130 and 24 detected events with a muon plus an up-going pion, respectively. 
The number of detected events Ndet(i,j, k) in each cell of Fig. H can be considered as 
the convolution (plus statistical fluctuations) of the reconstruction efficiency AE(i,j, k) 
with the unknown physical distribution of the events as a function of the variables in the 
parameter space V . To obtain these physical distributions, we evaluated the corrected 
number of events as the number of detected events in each cell divided by the recon- 
struction efficiency: N corr (i,j, k) = N det (i,j, k) / AE(i, j, k). Because of the uncertainties 
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Figure 8: a) Detector-unfolded distribution of the zenith scattering angle Aa of the pion with 
respect to the muon direction, for Aa > 45°. b) Detector- unfolded distribution of the up-going 
pion range inside the detector. Open points: pions with Aa > 45°; full points: pions with 
Aa > 60° (b). Only statistical errors are shown. The number of events corresponds to 1.55 y 
of data taking. 



in the reconstruction of Aa and from the measured track parameters, there is a 
small probability that an event in the (l,m,n) cell of (R w , Aa, D) will be measured in 
the (i,j,k) cell. This probability, which results in a smearing of the N det {i,j,k) distri- 
bution, has been evaluated using the simulated sample. The mean detection efficiency 
for Aa > 45°, averaged over the other parameters of V, is found to be 

t T = £ N det {i,j, k)/ ]T N corr (i,j, k) = (21 ± 2 stat ± A sys )% (4) 

i,j,k 

The systematic uncertainty is estimated by examining the effect of changing the sizes of 
the (R n , Aa, D) cells, and by evaluating the difference in the results obtained using the 
two hadronic interfaces to GEANT. 

From the distribution N corr {i, j,k), we obtain the detector-unfolded distribution of 
the pion-muon zenith scattering angle Aa and of the pion range inside the detector R n . 
The two distributions, which correspond to the x and y projection of J2k N CO rr(h j, k), 
respectively, are shown in Fig. [|. Given the low probability for events to be detected at 
small pion emission angles, the Aa distribution is shown for Aa > 45°. 

Fig. H shows the angular distribution of up-going charged particles versus cos($), 
where d is the zenith angle in the detector frame. The distribution is normalized to the 
total number of down-going single muons (12.2 x 10 6 ), and is corrected for reconstruction 
efficiency. 
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Figure 9: Angular distribution of the charged particles emerging from the floor, $, relative 
to the vertical direction. The flux is normalized to the total number of detected single muons 
(12.2 x 10 6 ). Open stars: uncorrected data. Full stars: corrected and detector-unfolded data 
for pions with scattering angles Aa > 45°. The pion energy threshold is around 100 MeV. 
Only statistical errors are shown. 



Using the detected number of muon plus up-going pion events and the mean detection 
efficiency quoted above, the estimated number of up-going pions emerging from the floor 
at the MACRO depth with angle of scattering Aa > 45° with respect to the muon 
direction and range > 25 g cm~ 2 , per down-going muon is found to be 

N n/lt = ^ = (io±l 5 ) icr 5 (5) 



6. A Model of High Energy Muon Interactions 

To compare with the results of the previous section, estimates of pion production by 
underground muons have been made using the hadronic interaction generator FLUKA[[8| 
combined with a model of hard muon scattering with the muon photonuclear cross 
sections of Ref. |J. The implications of this process for muon energy loss are discussed 
in ||10|| . The muon photonuclear cross section for the average muon energy at the MACRO 
depth is o" 7 (-E^ ~ 280 GeV) ~ 0.40 mb. This corresponds to a probability of 0.3% for 
a 300 GeV muon to have one hard scattering in one meter of rock. o 1 does not depend 
strongly on energy. For example, o^iE^ = 1000 GeV) ~ 0.55 mb. 

In || the interaction cross section for a muon of a fixed energy is given in terms 
of the fractional energy loss and of the q 2 transferred. Hadron production by muons 



is computed in the framework of the Weizsacker- Williams approximation [11 for the 
radiation of virtual photons, and the vector dominance model. Measured photon-vector 
meson couplings are used. Some uncertainties and limitations of the simulation arise 



from the fact that deep inelastic scattering is not included in FLUKA. However, this 
process is dominated by low q 2 interactions. In addition, although bremsstrahlung, pair 
production and other electromagnetic interactions of muons are considered, the transport 
of secondary e + e~ and 7 is not performed. Therefore, electromagnetic showers are not 
produced, nor is photo-production by real photons considered. 

The differential energy spectrum G(E^, h) of muons at the MACRO depth, along with 
the interaction model discussed above, have been used to estimate the flux of hadrons 
produced by muons at MACRO. An analytic approximation of G(E fJj ,h) |12[ is used, 
which assumes that the atmospheric inclusive muon spectrum above 1 TeV is described 
by the power law = KE' 1 , and gives a differential muon spectrum underground 

at the effective average depth h of G(E IM , h) = Ke~^'- l) (e^ + e (l - e^))" 7 ', with 
7=3.7, (3=0A (km w.e.)~ l and e=540 GeV jl2| . The average value of the rock overburden 
for MACRO, weighted by the measured muon flux, is h = 3.8 km.w.e.. Due to the shape 
of the rock overburden, there is a correlation between residual energy and angle. This 
correlation is neglected in the model. 

To simulate the production of up-going pions under MACRO, the muon energy is 
extracted from the energy spectrum shown above in the energy range 1 GeV < E < 
10 TeV. An interaction of the muon in the rock floor is then modeled, and the energy 
spectrum of the up-going particles exiting the floor is determined. All possible muon 
interactions are allowed, and the full hadronic cascade development is considered. For 
each particle type, the yield as a function of the kinetic energy and emission angle is 
then determined. We find that, as expected, the charged particle yield is dominated by 
charged pions, and that the differential pion spectrum decreases sharply with the kinetic 
energy and the emission angle Aa. The charged kaon contribution is found to be lower 
than the pion contribution by about one order of magnitude. Protons (and neutrons) 
have an energy spectrum which is considerably softer than that of pions. 

Next, in order to determine the number of muon plus up-going pion events due to 
this process observed in MACRO, muons with the energy distribution shown above are 
generated in (45 x 180) cells of (i^, (p^) using the intensity/depth relation given in [ |T3|j , 
and the slant depth obtained using a rock map of the Gran Sasso mountain (see ||13|| ). 
The measured density of the Gran Sasso rock (p = 2.657 g/cm?) is used in place of 
the usual standard rock. The absolute muon flux observed by MACRO provides the 
overall normalization for the simulation. Each simulated atmospheric muon undergoes 
a hard scattering in the rock around the apparatus. Only pion production is considered, 
according to the differential pion yield ^(E^, Aa) obtained by the simulation proce- 
dure above. The generated pions are then propagated to the apparatus. A simulated 
sample equivalent to 1.2 x 10 8 detected single muons was generated, corresponding to 
N y = 14.4 y of full detector livetime. 

Taking into account the real probability for hard scattering by an atmospheric muon, 
the predicted rate of muon plus up-going pion events in which both particles enter in 
the detector is (8 ±3) 10~ 5 ir/fi. The quoted error includes the uncertainty in the muon 
flux, the pion yield and the pion propagation in the detector. This result is in agreement 
with the measured value, quoted in Section 5. 



7. Evaluation of the Background due to Up-going 
Pions in Neutrino Studies. 



The flux of atmospheric muon neutrinos $ v (u) is often inferred from a measurement of 
the flux of up-going muons. MACRO can measure m the 100 GeV region using 

through-going muons and, in the 5 — 10 GeV region, using stopping muons. An up-going 
pion induced by an interaction of a down-going atmospheric muon in the rock around 
the detector results in a background event in the neutrino measurement if the muon is 
undetected. 

The rate of events in which the down-going muon is not detected, and in which the up- 
going pion simulates a neutrino-induced muon, has been evaluated using the simulation 
described in the previous section. In the simulated sample, equivalent to 14.4 y of data 
taking, 63 events are seen in which the up-going pion crosses two scintillator layers and 
the muon is undetected. Of these, 31 events contain pions which traverse more than 



200 g cm 2 of detector material. These events satisfy all of the criteria |14[ for an upward 
through-going muon. This gives a background rate due to up-going pions of 

Upward thoroughgoing [i background = (2.1 ± 0A stat ± 0.7 sys ) events /year 

The expected number of detected up-going muons from interactions in the rock be- 
low MACRO is approximately 200 events / y ear giving a contamination from this 
background of approximately 1%. The distribution of the cosine of the zenith angle of 
these background events is essentially flat. 

Partially contained muon neutrino interactions in the apparatus, giving rise to a 
down-going muon, and up-going stopping muons induced by a interaction in the rock 
below the apparatus, are identified by applying appropriate constraints to the observed 
tracks |TH]. In the following, these events are called stopping muons. The background 



for the stopping muon sample is evaluated by applying the analysis of ref. |TJj to the 
simulated events. 127 events are found in which the up-going pion satisfies the criteria 
used in this analysis for the identification of stopping muons, and in which the muon is 
not observed. This corresponds to a rate of these background events of 

Upgoing stopping fi background = (8.8 ± 0.9 sta t ± 0.8 SJ , S ) events / year 

The expected rate of stopping muons is about 80 events /year in MACRO, implying a 
contamination by up-going pions in this sample of about 10%. 

The background rates obtained above are specific to the MACRO detector. For 
smaller detectors, the probability of detecting the pion while missing the muon in- 
creases. For example, for a single MACRO supermodule the background/signal ratio for 
the through-going sample is almost twice the value quoted above. Other underground 
experiments measuring neutrino-induced muons should have a similar contamination 
due to up-going pions, with the exact background rate depending on the geometrical 
configuration of the detector, its depth, and on the momentum threshold applied to the 
up-going muons. To our knowledge, no experiment measuring the flux of atmospheric 



neutrinos using up-going muons (see references in |T4[]) has taken into account this back- 
ground. In particular, shallower experiments for similar energy threshold should have a 
significant contamination from these events, because the average muon energy, E^(h), 
and the corresponding pion yield, Y n /^, decreases slowly with decreasing depth, while 
the total muon flux increases exponentially. 



8. Conclusion 



In this paper we report the first measurement of up-going particles associated with down- 
going atmospheric muons in an underground detector. This measurement was possible 
because of the large area, good tracking and good time measurement capabilities of the 
MACRO apparatus. 

In 12.2 x 10 6 downgoing single muons we find a total of 243 events with an identifiable 
up-going particle spatially and temporally associated with a down-going muon. We 
have analyzed the experimental distributions of i) the total range of the pion from the 
production point to the stopping point, ii) the pion range inside the detector, Hi) the 
7r — \i zenith scattering angle and iv) the radial distance between the two particles in 
the bottom layer of scintillators. 

The probability A of events to be accepted, the tracking reconstruction efficiency E 
and the reconstruction efficiency AE in cells of the parameter space V(Rtt, Act, Aip, D) 
have been evaluated using a Monte Carlo simulation of the detector. This allows the cal- 
culation of the mean detection efficiency e r for events in which the 7r — \x scattering angle 
Aa > 45°. We find e r = (21 ± 2 stat ± A sys )%. We have presented the detector-unfolded 
distributions of the n — \x zenith scattering angle Aa, of the detector material thickness 
crossed by the upgoing particles, and of the angular distribution of charged particles 
emerging from the floor in the detector frame. The measured rate of up-going pions 
with range R n > 25 g cm" 2 per down-going muon at the Gran Sasso depth, corrected 
for the reconstruction efficiency of the MACRO apparatus, is N^u = (IOI2.5) 10~ 5 . 

A comparison of the data with a Monte Carlo calculation of up-going pion production 
in hard muon interactions shows good agreement in the measured pion production rate 
at the MACRO depth. 

The up-going particles generated by interactions of down-going muons constitute 
a background source for neutrino induced events. For MACRO, the background rate 
is (2.1 ± 0.4 stai ± 0.7 sys ) events / year for the upward through-going sample with a 
muon energy threshold of ~ 400 MeV, corresponding to 200 g cvrT 2 of detector. It is 
(8.8 ±0.9 stat ±0.8 SJ/S ) events I year for the up-going stopping sample. This background is 
of the order of 1% of the through-going muons and of the order of 10% of the up-going 
stopping muons. 
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